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We study switching between period-two states of an underdamped quantum oscillator modulated 
at nearly twice its natural frequency. For all temperatures and parameter values switching occurs via 
quantum activation; it is determined by diffusion over oscillator quasienergy, provided the relaxation 
rate exceeds the rate of interstate tunneling. The diffusion has quantum origin and accompanies 
relaxation to the stable state. We find the semiclassical distribution over quasienergy. For T — 0, 
where the system has detailed balance, this distribution differs from the distribution for T — » 0; 
the T = distribution is also destroyed by small dephasing of the oscillator. The characteristic 
quantum activation energy of switching displays a typical dependence on temperature and scaling 
behavior near the bifurcation point where period doubling occurs. 

PACS numbers: 03.65.Yz, 05.60.Gg, 05.70.Ln, 74.50.+r 



I. INTRODUCTION 

Switching between coexisting stable states underlies 
many phenomena in physics, from diffusion in soUds 
to protein folding. For classical systems in ther- 
mal equilibrium switching is often described by the 
activation law, with the switching probability being 
cx exp (— A?7/fcr), where AU is the activation en- 
ergy. As temperature is decreased, quantum fluctua- 
tions become more and more important, and below a 
certain crossover temperature switching occurs via tun- 
neling [J, [3, [l] . The behavior of systems away from ther- 
mal equilibrium is far more complicated. Still, for classi- 
cal systems switching is often described by an activation 
type law, with the temperature replaced by the charac- 
teristic intensity of the noise that leads to fluctuations 
i, i, i, 0, i, M, 0, EH, El- Quantum nonequilibrium 
systems can also switch via tunneling between classically 
accessible regions of their phase space [H, [H, [l^ [l^ . 

In addition to classical activation and quantum tun- 
neling, nonequilibrium systems have another somewhat 
counterintuitive mechanism of transitions between stable 
states. We call this mechanism quantum activation and 
study it in the present paper. It describes escape from 
a metastable state due to quantum fluctuations that ac- 
company relaxation of the system [l7|. These fluctua- 
tions lead to diffusion away from the metastable state 
and, ultimately, to transitions over the classical "bar- 
rier" , that is, the boundary of the basin of attraction to 
the metastable state. 

We study quantum activation for periodically modu- 
lated systems. Switching mechanisms for such systems 
are shown schematically in Fig. [U where g{Q) is the 
effective potential of the system in the rotating frame. 
This figure describes, in particular, a nonlinear oscil- 
lator studied in the present paper, which displays pe- 
riod doubling when its frequency is periodically modu- 
lated in time. The energy of periodically modulated sys- 
tems is not conserved. Instead they are characterized by 
quasienergy e. It gives the change of the wave function 
ipeit) when time is incremented by the modulation pe- 



riod Tp, ipeit + Tp) = eyip{—ieTF / h)ip^{t) , and is defined 
modulo 2'Kh/Tp fiir for. 

Coupling to a thermal reservoir leads to transitions be- 
tween the states of the system. For T = the transitions 
are accompanied by the creation of excitations in the 
thermal reservoir. The energy of the system decreases 
in each transition. However, the quasienergy may de- 
crease or increase, albeit with different probabilities VF| 
and VF|. This effect has quantum origin. It is due to the 
functions '0e being superpositions of the eigenfunctions 
|iV) of the energy operator of the system in the absence 
of modulation. Therefore bath-induced transitions down 
in energy |A^) — *■ — 1) lead to transitions V'e i^e' 




FIG. 1: (Color online) The effective double- well potential 
g[Q) of a parametrically modulated oscillator. Sketched are 
scaled period-two quasienergy levels (see Sec. |ll| in the ne- 
glect of interwell tunneling. The minima of g correspond to 
classically stable states of period two motion. The arrows in- 
dicate relaxation, diffusion over quasienergy away from the 
minima, and interwell tunneling. The effective Hamiltonian 
g{P,Q) is defined by Eq. ©, and g{Q) = g{P = 0,Q); the 
figure refers to ^ = in Eq. ((Gjl. 
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with e' ^ e. The values of W|, Vl^j^ for the latter transi- 
tions are determined by the appropriate overlap integrals 
and depend on e, e'. 

More probable transitions determine in which direc- 
tion, with respect to quasienergy, the system will most 
likely move. Such motion corresponds to relaxation over 
quasienergy. Fig. [prefers to the case Wi > W|. In 
this case relaxation corresponds to quasienergy decrease. 
The minima of g{Q) are the classical stable states. How- 
ever, quantum transitions in which quasienergy increases 
have a nonzero probability even for T = 0. Because 
this probability is less than Wi , such transitions lead to 
diffusion over quasienergy e away from the minima of 
g{Q). In turn, the diffusion leads to a finite-width dis- 
tribution over e and ultimately to activated-type over- 
barrier transitions between the wells in Fig. [T] In fact, 
discussed in this paper and sketched in Fig. [1] are period- 
two quasienergy states, with quasienergy e defined by the 
condition V'e(i + 2Ti?) = cxp{—2ieTp/h)il;g{t). They arc 
more convenient for the present problem; their relation 
to the standard quasienergy states is explained in Sec.HIl 

Of interest for the problem of switching is the semi- 
classical situation where the basins of attraction to the 
metastable states (the wells in Fig. [IJ have a large num- 
ber N of localized states. In this case the rate of tunnel- 
ing decay is exponentially small. The activation rate is 
also exponentially small since it is determined by the ra- 
tio of transition probabilities W| /Wi raised to the power 
N. Both the tunneling and activation exponents are oc N 
for the situation sketched in Fig. [T] Indeed, the tunnel- 
ing exponent is given by the action, in the units of h, for 
classical motion in the inverted effective potential —g{Q) 
from one maximum of ~~g{Q) to the other. It is easy to 
see that this action is of order N, unless g{Q) has a spe- 
cial form. Therefore for ^ 1 the activation exponent 
is cither much larger or much smaller than the tunneling 
exponent. 

In this paper we develop a theory of the statistical dis- 
tribution and consider switching of a parametrically mod- 
ulated quantum oscillator. We show that, irrespective of 
temperature, the activation exponent is smaller than the 
tunneling exponent. Therefore switching always occurs 
via activation, not tunneling, as long as the relaxation 
rate exceeds the tunneling rate. 

We study a nonlinear oscillator with frequency mod- 
ulated at nearly twice the natural frequency ujq- When 
the modulation is sufficiently strong, the oscillator has 
two stable vibrational states with periods 2tp w 271 /loq. 
These states differ in phase by tt but otherwise arc iden- 
tical [l^ . They correspond to the minima of the effective 
potential in Fig. [TJ The lowest quantized states in Fig.[T] 
are squeezed. Squeezing in a parametric oscillator has 
attracted interest in many areas, from quantum optics 
(l9l . [20l . [2l| to phonons microcantilevers [2^ and 

electrons and ions in Penning traps [l^l . Recent progress 
in systems based on J osep hson junctions and nanoelec- 
tromcchanical systems |25l. [26l[27l [28j makes it possible to 
study squeezed states in a well controlled and versatile 




-0.1 -0.05 



g 

FIG. 2: (Color online) The scaled inverse temperature R' = 
dR/dg of the distribution over scaled period-two quasienergy 
g of a parametrically modulated oscillator for different os- 
cillator Planck numbers n in the limit of a large number of 
intrawell quasienergy states. The figure refers to /i = —0.3. 



environment. Both classical and quantum fluctuations 
can be investigated and the nature of switching between 
the states can be explored. The results can be further 
used in quantum measurements for quantum computing, 
as in the case of switching between coexisting states of a 
resonantly driven oscillator p5| . 

The probability distribution and interstate transitions 
of a parametrically modulated oscillator have attracted 
considerable attention. Much theoretical work has been 
done for models where fluctuations satisfy the detailed 
balance condition either in the classical limit [2^ or for 
T = [33,[3l|,[33,[3l,[3i|. Generally this condition does 
not hold in systems away from equilibrium. In partic- 
ular a classical nonlinear parametric oscillator does not 
have detailed balance. Switching of such an oscillator 
was studied experimentally for electrons in Penning traps 
[35I . . The measured switching rate [11] agreed quan- 
titatively with the theory [13] ■ 

A quantum parametric oscillator also does not have 
detailed balance in the general case. The results pre- 
sented below show that breaking the special condition 
where detailed balance holds leads to a sharp change of 
the statistical distribution and the switching rate. This 
change occurs already for an infinitesimally small devia- 
tion from detailed balance, in the semiclassical limit. The 
fragility of the detailed balance solution is previewed in 
Fig. [2I This figure shows the effective inverse tempera- 
ture of the intrawell distribution over the scaled period- 
two quasienergy g. The function X^^R{g) is the expo- 
nent of the distribution X~^R(gn) = — In /9„, where p„ 
is the population of an n-th state and A oc 1/N is the 
scaled Planck constant defined in Eq. (jj]) below. The ef- 
fective inverse temperature dR/dg depends on g and 
differs from the inverse temperature of the bath T~^. The 
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FIG. 3: (Color online) Quantum activation energy of a tran- 
sition between the stable states of period two vibrations of a 
modulated oscillator (a phase-flip transition). The transition 
probability is Wa^f, oc exp(— _Ryi/A), where A, Eq. Q, is the 
effective Planck constant. The scaled switching exponent Ra 
is plotted as a function of the single scaled parameter /i that 
controls oscillator dynamics, Eq. for different values of 
the oscillator Planck number n. 



T = n = result for dR/dg is obtained from the solution 
with detailed balance (here n — [cx-p{hujo / kT) — 1]~ is 
the Planck number of the oscillator). It is seen from 
Fig. [2] that the T -> (n ^ 0) limit of the solution 
without detailed balance does not go over into the T = 
result. 

The effective quantum activation energy Ra « 
— AlogWsw, which gives the exponent of the switching 
rate Wsw, is shown in Fig. [31 It is equal to Ra = 
R{0) — R{gmin)i where g = and g = g^in are, respec- 
tively, the values of g{Q) at the top and the minimum of 
the wells in Fig. [1] The n = detailed balance value for 
Ra strongly differs from the fi —f value in a broad range 
of the control parameter fj, that characterizes the detun- 
ing of the modulation frequency from 2ujq, see Eq. (|5]) 
below. It is also seen from Fig. [3] that the quantum acti- 
vation energy Ra decreases with increasing temperature. 
Ultimately when the Planck number of the oscillator be- 
comes large, n ^ 1, we have Ra oc T^^, the law of ther- 
mal activation, and the results coincide with the results 
[sj obtained by a different method. 

In Sec. |TT]we obtain an effective Hamiltonian that de- 
scribes the oscillator dynamics in the rotating frame. We 
describe the classical dynamics and calculate semiclassi- 
cal matrix elements of the coordinate and momentum. In 
Sec. Inl] we derive the balance equation for the distribu- 
tion over quasienergy levels. This distribution is found in 
the eikonal approximation in Sec. II VI Its explicit form is 
obtained in several physically important limiting cases in 
Sec.|Vl In Sec. IVII the probability of interstate switching 



due to quantum activation is analyzed. In Sec. IVIII it is 
shown that the distribution obtained for T = 0, where 
the system has detailed balance, is fragile: in the limit of 
a large number of intrawell states it differs from the dis- 
tribution for T — > 0. It is also dramatically changed by 
even weak dephasing due to external noise. In Sec. IVIIII 
we discuss tunneling between the coexisting states of pe- 
riod two vibrations and show that interstate switching 
occurs via quantum activation rather than tunneling, if 
the relaxation rate exceeds the tunneling rate. In Sec. lIXI 
we summarize the model, the approximations, and the 
major results. 



II. DYNAMICS OF THE PARAMETRIC 
OSCILLATOR 

A. The Hamiltonian in the rotating frame 

We will study quantum fluctuations and interstate 
switching using an important model of a bistable sys- 
tem, a parametric oscillator. The Hamiltonian of the 
oscillator has a simple form 

Ho = + K + FcoaiiJFt)) + ^^q'' . (1) 

We will assume that the modulation frequency ojp is close 
to twice the frequency of small amplitude vibrations loq, 
and that the driving force F is not large, so that the 
oscillator nonlinearity remains small, 

5(jj = ^ujp - uJa , \Suj\<^uio, (2) 
|7|(g')«^o'- 

Here, (g^) is the mean squared oscillator displacement; 
in what follows for concreteness we set 7 > 0. 

We will change to the rotating frame using the 
canonical transformation U{t) = cxp {^—ia^doj pt/ 2^, 
where and a are the raising and lowering operators, 
d~ {fujjF)~^^'^ {ip + ujfi/^). It is convenient to intro- 
duce the dimensionless coordinate Q and momentum P, 

U\t)qU{t) = C[Pcos{ujFt/2)-Qiiiii{uJFt/2)] , (3) 
U\t)pU{t) = [Psm{LUFt/2) + Qcos{LUFt/2)] , 

where C = (2F/37)^/^. The commutation relation be- 
tween P and Q has the form 

[P, Q] = -iX, A = 3-/h/FujF ■ (4) 

The dimensionless parameter A will play the role of h 
in the quantum dynamics in the rotating frame. This 
dynamics is determined by the Hamiltonian 

Ho ^ U^HU - ihU^U w — g, (5) 

67 
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with 



.9 = 9{P,Q) 
1 



(l + /i)Q2. (6) 



Here we used the rotating wave approximation and dis- 
regarded fast oscillating terms oc exp(±mwj?t) , n > 1. 
As a result Hq is independent of time. 
In the "slow" dimensionless time 

T = tF/2ujF (7) 

the Schrodinger equation has the form iXdip/dr = gip. 
The effective Hamiltonian g, Eq. depends on one 
parameter 



fi = 2luf5uj/F . 



(8) 



For fi > —1, the function g{P, Q) has two minima. They 
are located at P = 0, Q = ±{n + 1)^/^, and gmin = 
— (/z + 1)^/4. For /i < 1 the minima are separated by a 
saddle at P = Q = 0, as shown in Fig. H) When friction 
is taken into account, the minima become stable states 
of the parametrically excited vibrations in the classical 
limit [ij, [we note that the function G in Ref. [13 is 
similar to g{P,Q), but has opposite sign]. 




FIG. 4: (Color online) The scaled effective Hamiltonian of 
the oscillator in the rotating frame g{P,Q), Eq. (|6]), for jj, = 
0.2. The minima of g{P, Q) correspond to the stable states of 
period two vibrations. 



The function g{P, Q) is symmetric, g{P, Q) = 
g{—P,—Q). This is a consequence of the time transla- 
tion symmetry. The sign change (P, Q) — > (— P, —Q) 
corresponds to the shift in time by the modulation pe- 
riod t t + 2tt/u!p, see Eq. ([3|). In contrast to the 
standard Hamiltonian of a nonrelativistic particle, the 
function ^(P, Q) docs not have the form of a sum of ki- 
netic and potential energies, and moreover, it is quartic 
in momentum P. This structure results from switching to 
the rotating frame. It leads to a significant modification 
of quantum dynamics, in particular tunneling, compared 
to the conventional case of a particle in a static potential. 



B. Quasienergy spectrum and classical motion 

1. Period-two quasienergy states 

The quasienergies of the parametric oscillator 
are determined by the eigenvalues gm of the operator 
g = g(P,Q), with an accuracy to small corrections 
^ F/ujp. They can be found from the Schrodinger equa- 
tion ^|77i)^°^ = gm\'m)^^\ Because of the symmetry of g 
the exact eigenfunctions \m)^^^ are either symmetric or 
antisymmetric in Q. 

The eigenstates |m)'^°-' give the Floquet states of the 
parametric oscillator ipeit + Tp) = cxjp{—ieTp / K)ip^{t) . 
From the form of the unitary operator U (t) = 
exp{-ia^au;pt/2) it is clear that U{Tp)\m)'-°'^ = ±|m)(°), 
with for symmetric and "— " for antisymmetric states 
|m)(°). Therefore we have e„i = [F"^ /&^)grn for symmet- 
ric and e„i ~ (F'^ /6^)gjn + {hujp/2) for antisymmetric 
states. 

It is convenient to define period-two quasienergy by 
the condition ipf{t + 2tp) = cxp{-'2ieTp / h)ip^{t) . In 
this case the relation e,„ = {F^/6^)gm holds for both 
symmetric and antisymmetric states. At the same time, 
there is simple one to one correspondence with the stan- 
dard quasienergy states. Period- two quasienergy states 
are particularly convenient for describing intrawell states 
\n) for small tunneling between the wells of g{P,Q). 
Such states are combinations of symmetric and antisym- 
metric states \m)^^\ For this reason we use period- two 
quasienergies throughout the paper. 

2. Classical intrawell motion 

We will study the intrawell states \n) in the semiclas- 
sical approximation. The classical equations of motion 



dQ ^ dP ^ 

^-dpg, ^ = -dQg 



(9) 



can be explicitly solved in terms of the Jacobi elliptic 
functions. The solution is given in Appendix A. There 
arc two classical trajectories for each g < 0, one per well 
of g{P,Q) in Fig. [51 They are inversely symmetrical in 
phase space and double periodic in time. 



Q(T + rp) = 0(r), P{T + Tp)^PiT). 



(10) 



with one real period Tp^^ and one complex period Tp^' 



This means that t„ 



0,±1, 



and 



(1) 



n2Tp'^'^ , where ni. 



rW^2V2|5|-i/4i^(^,), 
f)=z2V2|g|-i/4^'(^,). 



(11) 



Here K{mj) is the complete elliptic integral of the first 
kind, K'{mj) = K[l — mj), and mj = mj{g), 

, , (M+l-2|g|^/^) (/.-l + 2|.g|V^) 
mj(<?) = *0 (12) 
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is the parameter [38| . The real part of m j can be positive 



has not only 



or negative; for Remj < 0, the period Tp 
imaginary, but also a nonzero real part. 

The vibration frequency is = 27r/rp 

2Tr/Tp^\g). It monotonically decreases with increasing g 
in the range where g{P, Q) has two wells, — (/x + 1)^/4 < 
g < 0, and Lo{g) for g 0, i.e. for g approaching 
the saddle-point value. The periodicity of Q{t), P{t) in 
imaginary time turns out to be instrumental in calculat- 
ing matrix elements of the operators P, Q on the intrawell 
wave functions In) as discussed in the next section. 



The semiclassical matrix elements 



Im ( 



' 1 

m< 
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Of central interest to us will be the case where the 
number of states with gn < 0, that is the number of 
states inside the wells of g{P,Q) in Fig.[4l is large. For- 
mally, this case corresponds to the limit of a small effec- 
tive Planck constant A. For A ^ 1 the wave function |n) 
in the classically accessible region can be written in the 
form 



|n) =C'(ap.9)-i/2cxp [tSn{Q)/X]+c.c., 
rQ 

SniQ) - / P{Q',gn)dQ' . 



(13) 



Here, C" is the normalization constant, P{Q,gn) is the 
momentum for a given Q as determined by the equa- 
tion g(P,Q) = gn, Qn is the classical turning point, 
P{Qn,gn) = 0, and the derivative dpg is calculated for 
P = PiQ,gn). 

We parameterize the classical trajectories Q by their 
phase (j) = uj{g)T, which we count (modulo 2tt) from its 
value at Qn- With this parametrization we have 



(Q) - SniQ) ^ Xmcj) , 



(14) 



for \m\ <ti A~^; on a classical trajectory, Q is an even and 
P is an odd function of (f>. 

With Eqs. and (|14p . we can write the semi- 

classical matrix element of the lowering operator a ^ 
(2A)~^/^(P — iQ) as a Fourier component of a((j);g) = 
(2A)~^/^[P(r; q) — iQ{T;g)] over the phase = i^{g)T 



{n + TO|a|n) 



(15) 



— / d(j)e^p{-im(f>)a{(l);g) , 



to lowest order in A. As we discuss below, these matrix 
elements determine relaxation of the oscillator. 

We will be interested in calculating am{g) for 5 < 0. 
Then if we neglect interwell tunneling, we have two sets of 
wave functions |n), one for each well. The interwell ma- 
trix elements of the operator a are exponentially small. 
The intrawell matrix elements are the same for the both 



FIG. 5: (Color online) The contours of integration in the 
complex phase plane for calculating Fourier components of the 
oscillator coordinate Q and momentum P. The left and right 
panels refer, respectively, to the negative and positive real 
part of the parameter of the elliptic functions mj, Eq. (|12p 
[the imaginary part of mj is infinitesimally small]. Both Q 
and P have poles at (/> = 4>*,4'**i but a(0; g) has a pole only for 
4> = (j),. The increment of phase by (f>o leads to a transition 
between the trajectories with a given g in different wells of 



wells except for the overall sign, because Q and P in dif- 
ferent wells have opposite signs on the trajectories with 
the same g. We will consider the matrix elements am{gn) 
for the states in the right well Q > 0. 

The integral ([TS]) can be evaluated using the double 
periodicity of the functions Q, P. It is clear that a((/) -I- 
2-K, g) = a((j)', g) for any complex (f>. In addition, from 
Eq. (|A3P for g < we have a{(j) + 0o,5) = ^a(0;5), 
where (f>Q — 7r(l -I- t^"^ /t^'^) is a complex number, which 
is determined by the ratio of the periods of Q and P. 
The shift in time by 4>Qlu}{g) corresponds to a transition 
from Q{t),P{t) in one well to Q{t),P{t) in the other 
well. We can now write the matrix element am{g) as 



am(g) 



1 

2^ 



[1 



d0e-™^a(0;g), (16) 



where integration is done along the contour C in Fig. [51 
It is explained in Appendix |^ that each of the functions 
Q ((l)/Lu(q)), P {4)/oj{g)) has two poles inside the contour 
C [331 . Therefore it is easy to calculate the contour inte- 
gral The result has a simple form 



a^{g) = -^(2A)-l/2^(.9) 



cxp(— im0,) 
1 + exp(— im(/)o) 



with given by the equation 



Imc 



<-Im 



l + M-2|.g|i/^ 



■^(2)/^(l) 



(17) 



(18) 



-Im ( 
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An important feature of the matrix elements am{g) is 
their exponential decay for large |77i|. From Eqs. (|17p and 

m 

o-mig) oc exp [— mlni(0o - 4>*)] , m » 1 , (19) 
am(.9) oc exp[— |m|lm0*] , — m ^ 1 . 

We note that the decay is asymetric with respect to the 
sign of m. This leads to important features of the prob- 
ability distribution of the oscillator. 

III. BALANCE EQUATION 

Coupling of the oscillator to a thermal reservoir leads 
to its relaxation. We will first consider the simplest type 
of relaxation. It arises from coupling linear in the oscil- 
lator coordinate q and corresponds to decay processes in 
which the oscillator makes transitions between neighbor- 
ing energy levels, with energy « hujo transferred to or 
absorbed from the reservoir. We will assume that the os- 
cillator nonlinearity is not strong and that the detuning 
duj of the modulation frequency is small, whereas the den- 
sity of states of the reservoir weighted with interaction is 
smooth near loq- Then the quantum kinetic equation for 
the oscillator density matrix p in the rotating frame has 
the form 

^=^X~^[p,g]-f,p, (20) 

VP = V [{tt- + l)(a^ap — 2apa) + pa) a) 
+n{aa) p — 2a) pa + paa))\ , 

where rj is the dimensionless relaxation constant and n = 
[exp (huJo/kT) — is the Planck number. 

We will assume that relaxation is slow so that the 
broadening of quasienergy levels is much smaller than 
the distance between them, rj ^{g)- Then off-diagonal 
matrix elements of p on the states |n) are small. We note 
that, at the same time, off-diagonal matrix elements of 
p on the Fock states of the oscillator \N) do not have to 
be small. 

To the lowest order in rj/uj{g) relaxation of the diagonal 
matrix elements pn = {n\p\n), is described by the balance 
equation 

^ = -277^ [Wnn'Pn " W„.„p„0 , (21) 

n' 

with dimensionless transition probabilities 

Wnn' = (n+l)|(n'|a|7i)p+,l|(n|a|n')p. (22) 

It follows from Eqs. HT]) and ([22]) that, even for T = 0, 
the oscillator can make transitions to states with both 
higher and lower gr, with probabilities Wnn' with n' > n 
and n' < n, respectively. This is in spite the fact that 
transitions between the Fock states of the oscillator for 
T = are only of the type \N) ^\N -1). 



The explicit expression for the matrix elements (fT7|) 
makes it possible to show that the probability of a tran- 
sition to a lower level of g is larger than the probability 
of a transition to a higher level, that is, Wn'n > Wnn' 
for n' > n. Therefore the oscillator is more likely to 
move down to the bottom of the initially occupied well 
in Fig. m This agrees with the classical limit in which 
the stable states of an underdamped parametric oscillator 
arc at the minima of g{P,Q). However, along with the 
drift down in the scaled period-two quasienergy g^ even 
for T = there is also diffusion over quasienergy away 
from the minima of g{P,Q), due to nonzero transition 
probabilities Wnn' with n' > n. 

Equation ([2T|) must be slightly modified in the presence 
of tunneling between the states with gn < 0. The mod- 
ification is standard. One has to take into account that 
the matrix elements of p depend not only on the number 
n of the period-two quasienergy level inside the well, but 
also on the index a which takes on two values a = ±1, 
that specify the wells of g{P, Q). These values can be as- 
sociated with the eigenvalues of the pseudospin operator 
Uz- The matrix elements of the operator f]p on the wave 
functions of different wells are exponentially small and 
can be disregarded. The operator describing interwell 
tunneling can be written in the pseudospin representa- 
tion as —i{2X)~^T[g) [ux, p], where T{g) is the timncling 
splitting of the states in different wells. 

We assume that the tunneling splitting is much smaller 
than rj. Then after a transient time r ~ t?"^, there is 
formed a quasistationary distribution over the states |n) 
inside each of the wells of Q) in Fig.[31 This distribu- 
tion can be found from Eq. (|2ip using the wave function 
\n) calculated in the neglect of tunneling. Interwell tran- 
sitions occur over much longer time. 

IV. DISTRIBUTION OVER INTRAWELL 
STATES 

The stationary intrawell probability distribution can 
be easily found from Eqs. (pij) and if the number of 
levels with g„ < is small. Much more interesting is the 
situation where this number is large. It corresponds to 
the limit of small A. We will be interested primarily in the 
quasistationary distribution over the states in the well in 
which the system was initially prepared. It is formed over 
time ~ ri~^ and is determined by setting the right hand 
side of Eq. pT|) equal to zero. The resulting equation 
describes also the stationary distribution in both wells, 
which is formed over a much longer time given by the 
reciprocal rate of interwell transitions. 

For A ^ 1 we can use the Wcntzel-Kramers-Brillouin 
(WKB) approximation both to calculate the matrix el- 
ements in Wnn' ([22]) and to solve the balance equation. 
The solution should be sought in the eikonal form 

Pn = exp [-i?„/A] , Rn = R{gn)- (23) 
It follows from Eq. p9p that the transition probabilities 
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Wn+rnn rapidly decay for large 
balance equation (pij) we can set 



Therefore in the 



Pn+m ~ Pn CXp | 

R'{g) = dR/dg. 



-muj{gn)R'{gn)] 



(24) 



The corrections to the exponent in Eq. (|24p of order 
Xm^Lu'^R", \m?uju)'R' are small for A ^ 1 (prime in- 
dicates diflFerentiation over g). We note that iv{g)R'{g) is 
not small in the quantum regime and the exponential in 
p4|) will not be expanded in a series in R'. 

From Eq. pT|) . the function R'{g) is determined by the 
polynomial equation 



^w„+„,„(i-C) = 0, 

rn 

( = cxp [-u;{gn)R'ign)] ■ 



(25) 



Here the sum goes over positive and negative m. In ob- 
taining Eq. (j25p we have used the relation Wn+rnn — 



n n—m 



for 



< n, 



which is the consequence 



of the WKB approximation for the matrix elements: 
Aw(g)dln [a,„(.g)] /dg < 1 for A < 1. 

Equation ([25]) has a trivial solution ^ = 1, which is 
unphysical. Because all coefficients Wn+mn are positive 
and Wn+mn > Wn-rnn for m > 0, ouc Can show that 
the polynomial in the left hand side of Eq. (|25p has one 
extremum in the interval < ^ < 1. Since its derivative 
for ^ = 1 is negative and it goes to — oo for ^ — > 0, it has 
one root in this interval. This root gives the value i?' ((?„). 
Since ^ < 1, we have R'{g) > 0. In turn R'{g) gives 
R{g) and thus the distribution p„. In obtaining R{g) 
from R'{g), in the spirit of the eikonal approximation, 
one should set R{gmin) = 0, where gmin = —(/i +1)^/4 is 
the minimal quasienergy. 

It is seen from Eq. that X^^R'{g) has the mean- 
ing of the effective inverse temperature of the distribution 



-0.12 -0.06 -0.4 -0.2 



FIG. 6; (Color online) The scaled inverse temperature R' of 
the distribution over scaled period-two quasienergy g for two 
values of the control parameter fi and for different oscillator 
Planck numbers n. For n > 0.1 the function R' only weakly 
depends on g inside the well, gmin < g < 0. 



over period- two quasienergy. The function R'{g) can be 
calculated for different values of the control parameter p, 
and different oscillator Planck numbers n by solving Eq. 
(|25p numerically. The numerical calculation is simpli- 
fied by the exponential decay of the coefficients Wn+mn 
with |m|. The results are shown in Fig. [6l The fimc- 
tion R'{g)/\ smoothly varies with g in the whole range 
<7min < <7 < of intrawell values of g, except for very 
small n. 

Numerical results on the logarithm of the distribution 
R{gn) = — Alnpn obtained from Eq. (P5|) for n = 0.1 are 
compared in FiglTja) with the results of the full numer- 
ical solution of the balance equation ([21]) . In this latter 
calculation we did not use the WKB approximation to 
find the transition probabilities Wn n' ■ Instead they were 
obtained by solving numerically the Schrodinger equation 
g\n) = (?„|rt) and by calculating Wnn' as the appropri- 
ately weighted matrix elements (n|a|n') , Eq. ([22]) . 



In the calculation we took into account that, for small 
A, the levels oi g form tunnel-split doublets. The tunnel- 
ing splitting is small compared to the distance between 
intrawell states Xuj{g)- The exact eigenfunctions of the 
operator g are well approximated by symmetric and anti- 
symmetric combinations of the intrawell wave functions. 
This allows one to restore the intrawell wave functions 
from the full numerical solution and to calculate the ma- 
trix elements Wn n' ■ By construction, such numerical ap- 
proach gives R{g) only for the values of _g that correspond 
to quasienergy levels It is seen from Fig. [Tlja) that 
the two methods give extremely close results for small A 
(see also below). 
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FIG. 7: (Color online) Comparison of the results of the eikonal 
approximation for the scaled logarithm of the probability dis- 
tribution R over period-two quasienergy g (solid lines) with 
the results obtained by direct calculation of the transition 
probabilities followed by numerical solution of the balance 
equation (squares). 
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V. THE DISTRIBUTION IN LIMITING CASES 

The effective inverse temperature R'{g)/\ and tlie dis- 
tribution pn can be found in several limiting cases. We 
start with the vicinity of the bottom of the wells of 
g{P,Q). Here classical vibrations of P, Q are nearly 
harmonic. To leading order in dg = g — 5min we have 
{n + m|a|n) cx S\m\,i<>g^^^ for m ^ 0, and therefore the 
transition probabilities Wn+mn « S\„i^ i. Then Eq. (pS]) 
becomes a quadratic equation for ^, giving 

i?'(5mi„) = + (26) 

^ in + 2){2n + 1) + 2V7rn" 

The inverse effective temperature -R'((?min) as given by 
Eq. (pS)) monotonically decreases with the increasing 
Planck number n, i.e., with increasing bath temperature 
T. We note that i?'(5min) smoothly varies with n for 
low temperatures n <C 1, except for small Wc have 
dR'{g,nin)/dn = -~2^/^{fi + 1)/^ for n 0. For ^ 0, 
on the other hand, we have i?'(.gmin) ~ hujp/AkT, i.e. 
the effective inverse temperature for g = gmin is simply 
cx T-\ 

We now consider the parameter ranges where R'{g) can 
be found for all g. 

A. Classical limit 

For n 3> 1 the transition probabilities become nearly 
symmetric, \Wn+mn - Wn^mn\ < Wn+mn- As a result, 
the effective inverse temperature i?'/A becomes small, 
and the ratio pn+m/pn can be expanded in i?'(f/„)/A. 
This gives 

=2L^-i((?)^mW„+„„/^m2t^„+„„. (27) 

It is shown in Appendix [B] that Eq. can be written 
in the simple form 

^'(5) = ^^M(9)/Nig) , (28) 
M{g)^ ff dQdP, 

N{g) = \ll dQdP{dlg + dig) , 

where the integration is performed over the area A{g) of 
the phase plane (Q,P) encircled by the classical trajectory 
(5(r), P{t) with given g. 

From Eq. ((28|) . the effective inverse temperature i?'/A 
is oc [A(2n + 1)] . In the high-temperature limit {2n -t- 
1) « 2kT/hL0Q. Therefore R'/\ is oc and does 

not contain /i, as expected. Eq. (|28|) in this limit co- 
incides with the expression for the distribution obtained 
in Rcf. [13 using a completely different method. 



B. Vicinity of the bifurcation point 

The function R'{g) may be expected to have a simple 
form for p close to the bifurcation value = — 1 where 
the two stable states of the oscillator merge together and 
.g(P, Q) becomes single well. This is a consequence of 
the universality that characterizes the dynamics near bi- 
furcation points and is related to the slowing down of 
motion and the onset of a soft mode. The situation we 
are considering here is not the standard situation of the 
classical theory where the problem is reduced to fluctua- 
tions of the soft mode. The oscillator is not too close to 
the bifurcation point, its motion is not overdampcd and 
the intcrlevcl distance exceeds the level broadening. Still 
as we show R{g) has a simple form. 

The motion slowing down for p approaching —1 leads 
to the decrease of the vibration frequencies. From 

Eq. dni), 

u{g) ^2n/T^^W\g\'/^ <[{p + l) /2f" 

for /i + 1 <C 1. Therefore one may expect that Lo{g)R'{g) 
becomes small, and we can again expand Pn+m/Pn in 
uj{gn)R'{gn), as in Eq. l|7f)). One can justify this ex- 
pansion more formally by noticing that the transition 
probabilities Wn+mn are nearly symmetric near the bi- 
furcation point \(Wn+mn -Wn-mn)\ < Wn^mn- In 

the present case the latter inequality is a consequence 
of the relation Im ((/)o — 2(/)h.) ^ 1 in Eq. p?)) . which 
leads to ||a„m(5)| - |am(g)| |< \am{g)\- In turn, the 
above relation between 0o and 0, can be obtained from 
Eqs. (fT^ and which show that the right-hand side 
of Eq. (UHl) is « - |mj/ (1 - mj)|^/^ for /x-hK 1. Since 
|cn (A' + i X')! = \fnjl (1 — toj)|^^^ [sll, we have from 

Eqs. dni) and ^ Im0, w Imr^^V27-i^^ = Im</>o/2. 

It follows from the above arguments that near the bi- 
furcation point R'{g) is given by Eq. (|28p for arbitrary 
Planck number n, i.e., for arbitrary temperature. Since 
Sgf/ -I- dig « 2 for -I- 1 ^ 1, we obtain from Eq. (^5)) a 
simple explicit expression 

i?'(.9)«2/(2n+l), (29) 

Eq. (pg)) agrees with Eq. (pS)) near t^min in the hmit p. + 
1^1. It shows that the inverse effective temperature 
R' / A is independent of g. It monotonically decreases with 
increasing temperature T. 

C. Zero temperature: detailed balance 

The function R'{g) can be also obtained for T = n = Q. 
This is a consequence of detailed balance that emerges in 
this case [S^]. The detailed balance condition is usu- 
ally a consequence of time reversibility, which does not 
characterize the dynamics of a periodically modulated 
oscillator. So, in the present case detailed balance comes 
from a special relation between the parameters for T = 0. 
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Detailed balance means that transitions back and forth 
between any two states are balanced. It is met if the ra- 
tio of the probabilities of direct transitions between two 
states is equal to the ratio of probabilities of transitions 
via an intermediate state 

One can sec from Eqs. (|17l) and ([22]) that the condi- 
tion ((50]) is indeed met for n = 0. Therefore the balance 
equation has a solution Pn/ Pn' = Wn'n/Wnn' , which im- 
mediately gives 

i?'(g) = 2L^-i(5)Im(0o-20,). (31) 

The function R'{g) is plotted in Fig. [5] for two values 
of the control parameter fi. It displays different behavior 
depending on the sign of fi. For /i < the parameter 
of the elliptic function mj, Eq. p^ . is negative for all 

(1 2) 

g < 0. Therefore the periods Tp ' {g) of P{t), Q{t), 
Eq. ((TT|) . are smooth functions of g, except near jmin = 
where Imri'^ « (i + ^)-i/2 |in (g _ g^i„)|. 

This divergence of Imr^ is seen in Fig.[8l The function 
R'{g) remains finite near f/min, Eq- (ES])- 



15n . , ^ 1 15p ' , 

n = -0.3 ^ = 0.3 




01. . , . 1 QL, : , . 1 

-0.12 -0.06 -0.4 -0.2 

g g 

FIG. 8: (Color online) The scaled inverse temperature R' 

(2) 

and the imaginary part of the period Tp of the oscillator 
coordinate and momentum Q{t; g), P{t; g) for T = 0. For 
p> both R' and Imr^^^ have a logarithmic singularity. 

For p > 0, on the other hand, the function R'{g) has a 
singularity. Its location gd is determined by the condition 
mj = 0, which gives gd = —(1 — pY' /A. For small \g — 

(2) 

gd\ we have Imi/iQ oc ImTp cx |ln (jg — grf|)|, whereas 
(/)*, Eq. remains finite for g = gd- Therefore R'{g) 
diverges logarithmically at gd, as seen from Fig.[5{b). 

Physically the divergence of the inverse temperature 
is related to the structure of the transition probabili- 
ties Wnn' for n ~ 0. For g gd we have Wnn+m 
{gn - gdf"" for TO > 0, see Eqs. (HT]) and (EH). This 
means that, for gn close to gd, diffusion towards larger 



g slows down. The slowing down leads to a logarithmic 
singularity of R'{g). The eikonal approximation is inap- 
plicable for g close to gd- However, the width of the range 
of g where this happens is small, Sg ^ A, as follows from 
the discussion below Eq. ((24|) . In addition, there arc cor- 
rections to the balance equation due to off-diagonal terms 
in the full kinetic equation. These corrections give extra 
terms cx rf /uP'{gd) in the transition probabilities Wnn'- 
The analysis of these corrections as well as features of 
R that are not described by the eikonal approximation 
is beyond the scope of this paper because, as we show, 
these features are fragile. 

The function R{gn) = In pn obtained by integrating 
Eq. PT|) is compared in Fig. [7^b) with the result of the 
numerical solution of the balance equation (j2ip with nu- 
merically calculated transition probabilities Wn n' ■ The 
semiclassical and numerical results are in excellent agree- 
ment. We checked that the agreement persists for differ- 
ent values of the control parameter p and for almost all 
A <C 1 , except for a few extremely narrow resonant bands 
of A. 



VI. SWITCHING EXPONENT 

Quantum diffusion over quasienergy described by 
Eq. pT|) leads to switching between the classically stable 
states of the oscillator at the minima of g{P, Q) in Fig. H) 
The switching rate W^sw is determined by the probability 
to reach the top of the barrier of g{P, Q), that is by the 
distribution p.„ for such n that (7„ = 0. To logarithmic 
accuracy 

.0 

Wsw = QwXexp(-i?A/A), Ra= R'{g)dg, {32) 

where R'{g) is given by Eq. ([25]) . The parameter Csw is 
of the order of the relaxation rate rj due to coupling to a 
thermal bath. 

The quantity Ra plays the role of the activation energy 
of escape. The activation is due to quantum fluctuations 
that accompany relaxation of the oscillator, and we call it 
quantum activation energy. As we show in Sec. IVIIIi Ra 
is smaller than the tunneling exponent for tunneling be- 
tween the minima of g{P, Q). Therefore if the relaxation 
rate ry exceeds the tunneling rate, switching between the 
stable states occurs via quantum activation. 

Quantum activation energy Ra obtained by solving 
Eq. ([25]) numerically for R' is plotted in Fig. [3l It de- 
pends on the control parameter p and the Planck number 
fi, and it monotonically increases with increasing p and 
decreasing n. Close to the bifurcation point pB = — 1, 
i.e., for p — pb ^ 1, it displays scaling behavior with 
p- PB- From Eq. ^ 

RA = \{2n + l)-^{p-pBf, C = 2. (33) 

The scaling exponent C, — 2 coincides with the scaling 
exponent near the bifurcation point of a classical para- 
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FIG. 9: (Color online) Quantum activation energy of switch- 
ing between the states of parametrically excited vibrations 
for different oscillator Planck numbers n as a function of 
the scaled frequency detuning [i. The transition rate is 
Wsw oc exp(— i?,4/A). With increasing n, the value of Ka 
multiplied by 2?! + 1 quickly approaches the classical limit 
n ^ 1 shown by the dashed line. In this limit the ratio -Ra/A 
is oc and does not contain h. 



for T = docs not generally apply even for infinitesi- 
mally small but nonzero temperature, in the semiclassi- 
cal limit. This is a consequence of this solution being of 
singular nature, in some sense. A periodically modulated 
oscillator should not have detailed balance, because the 
underlying time reversibility is broken; the detailed bal- 
ance condition pop is satisfied just for one value of n and 
when other relaxation mechanisms arc disregarded. 

Formally, in a broad range of [i the correction oc n 
to the T = solution diverges. The divergence can be 
seen from Eqs. HH), (gl]), The transition proba- 

bilities Wnn+m havc tcrms cx h which vary with m as 
exp [— 2m Im(/)*] for m ^ 1. At the same time, the 
T = solution §^ gives ^'"^ = exp[muj{g)R'(g)] = 
exp [2 m lm(0o ^ 2(f>*)]- Therefore for the series (|25p with 
the term c>c n in Wn n+m to converge we have to have 



lm(0o - 30*) < 0. (35) 



metric oscillator where the oscillator motion is still un- 
derdamped in the rotating frame, i.e., /i — fiB is not to 
small [33|. It can be seen from the results of Rcfs. l40l [4ll 
that the exponent ( — 2 also describes scaling of the 
activation energy of escape due to classical fluctuations 
closer to the pitchfork bifurcation point where the motion 
is necessarily overdamped. 

In the classical limit 2n -I- 1 ^ 1 we have from Eq. (pS]) 
Ra oc (2n + 1) oc T^^, i.e., the switching rate Wsw oc 
exp [— (RAksT / X) /ksT], with temperature independent 
RAksT/X being the standard activation energy. The 
quantity i?yi(2n. + 1) in the classical limit as obtained 
from Eq. is shown with the dashed line in Fig. O 

It is seen that Ra quickly approaches the classical limit 
with increasing Planck number fi, so that even for n ~ 0.1 
the difference between RA{2n + 1) and its classical limit 
is < 15%. Effectively it means that, for small n 1, the 
exponent in the probability of quantum activation can 
be approximated by the classical exponent for activated 
switching in which one should replace 

ksT fujo ■ (34) 



The condition ([55)1 is met at the bottom of the wells 
of g(P, Q) and also close to the bifurcation values of 
the control parameter, — /i^ ^ 1. We found that, 
with increasing /x, the condition (j35p is broken first for 
g approaching the barrier top, 5 — > 0. In this region 

Im c/jojlm 0* cx^ |ln|5|| ■ A somewhat tedious calcu- 
lation based on the properties of the elliptic functions 
shows that the condition ([55)1 is violated when fi > —1/2, 
and for ^ = — 1/2 we have Im ((/jq — 30*) — > for g ^ 0. 
The increase in /i leads to an increase of the range of 
g where Eq. ([55)) docs not apply. The detailed balance 
distribution is inapplicable for T — > in this range. 

To find the distribution for small Planck number n in 
the range where Im (00 — 30*) > we seek the solution 
of the balance equation ([25)1 in the form 



R'ig) = 2c^-i(.g)[Im0*(.g) - e] , e « 1 . (36) 



VII. FRAGILITY OF THE DETAILED 
BALANCE SOLUTION 

It turns out that the expression for the distribution 
and (pij) found from the detailed balance condition 



This solution docs not give diverging terms for n > 0. 
The terms Wnn+mQy^"p[mu}{g)R' {g)] in Eq. (^5]) are oc 
nexp(— 2me) for m ^ 1, and their sum over to is oc n/e. 
Using the explicit expression for the transition probabil- 
ities we obtain that, to leading order in n 
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cxp [m lni(0o — 3(/>*)] sinh (m Imc/)*) / |cos (to Im0o/2)|^ 



(37) 



Here again, the sum goes over positive and negative to. In 
obtaining Eq. ([57]) we used the general expression (fT7|l for 
the coefBcicnts am{g) in Wn+mn- In fact, the asymptotic 
expressions (|19[) for a^ig) allow one to calculate the sum 
over TO in Eq. (j37p explicitly and the result provides a 
very good approximation for e. 

Equations ((36)) and ((37)) give the effective inverse tem- 
perature i?'/A and therefore the distribution p„ as a 
whole for fi <C 1. It is clear that the fi ^ limit is 
completely different from the detailed balance solution 
((3T)) for n = 0, that is, the transition to the T = n = 
regime is nonanalytic. The difference between the 7I, = 
and the n ^ solution is seen in Figs. [D and O 

It follows from Eq. ((37)) that the perturbation theory 
diverges for the value of g where Im {<f>Q — 30*) = 0. At 
such g the n — > solution ((55)) coincides with the n = 
solution ((nj) . For smaller g the condition ((55)) is satisfied 
and the T = solution for R'{g) applies; the corrections 
to this solution are oc n ^ 1 . The derivative of the effec- 
tive inverse temperature R'{g)/X over g is discontinuous 
at the crossover between n — > and n = solutions, 
as seen in Fig. [2] Figure [2] illustrates also the smearing 
of the singularity of R'{g) due to terms ex n in Wnn', 
which is described by the numerical solution of the bal- 
ance equation ()25p . Away from the crossover the analyt- 
ical solutions provide a good approximation to numerical 
results. 



A. Breaking of detailed balance solution by 
dephasing 

Dephasing plays an important role in the dynamics of 
quantum systems. It comes from fluctuations of the tran- 
sition frequency due to external noise or to coupling to 
a thermal reservoir. A simple mechanism is quasielastic 
scattering of excitations of the reservoir off the quantum 
system. Since the scattering amplitude depends on the 
state of the system, the scattering leads to diffusion of 
the phase difference of different states. 

For an oscillator, dephasing has been carefully stud- 
ied, both microscopically and phenomenologically, see 
Refs. HH, 0, and papers cited therein. It leads to 
an extra term ^if^p in the quantum kinetic equation 
(EO)), with 



(38) 



where 77^*^ is the dimensionless dephasing rate. 

If both ryph and 77 are small compared to oj{g)^ popu- 
lations of the steady states p„ are described by the bal- 
ance equation (|2ip in which one should replace rjWn'n 



(Ph) 
n'n ■ 



with 

(ph) _ 



(n a' a\n 



(39) 

n\. we have 



It follows from Eq. ()19|) that, for large | 

^r^l% exp [—2|n' — n|Im that is, the transition 
probability exponentially decays with increasing |n' — n|, 
and the exponent is determined by Im(/)*. 

Even slow dephasing is sufficient for making the de- 
tailed balance condition inapplicable. Mathematically, 
the effect of slow dephasing is similar to the effect 
of nonzero temperature. If the condition ()35p is vi- 
olated, the sum ^ W,j'^^j„ exp (— TOw(,g„)i?'(g„)) with 
R'{g) given by the detailed balance solution ()3ip diverges. 
The correct distribution for the appropriate g and fi is 
given by Eq. ((55)) . The parameter e = e{g) is given by 
Eq. ((57)) in which n is replaced, 

2 



exp (2iTO(/)*) / [exp [imcj)!)) + 1] 



(40) 



It is seen from Eq. ((10)) and Figs. [5] and [5] that, for low 
temperatures, even weak dephasing, rjP^/rj ^ 1, leads to 
a very strong change of the probability distribution. 

The fragility of the detailed balance solution discussed 
in this section is a semiclassical effect. It occurs if the 
number of states in each well TV oc is large. Formally 
we need eexp[2cA^Im (0o ~ •^4'*)] ^ li with c = c(g) ~ 1. 
In other words, the detailed balance solution is fragile 
provided A is sufficiently small. A full numerical solution 
of the balance equation confirms that, when A is no longer 
small, the distribution for small n, rj^^/r] is close to that 
described by the n = 0, rj^^/i] = 0. 



VIII. TUNNELING 

The oscillator localized initially in one of the wells of 
the effective Hamiltonian g{P, Q) can switch to another 
well via tunneling. For small A tunneling can be de- 
scribed in the WKB approximation. We will first find 
the tunneling exponent assuming that the oscillator is 
in the lowest intrawell state and show that it exceeds 
the quantum activation exponent Ra/X- We will then 
use standard arguments to show that this is true inde- 
pendent of the initially occupied intrawell state, for all 
temperatures. 

If we now go back to the original problem of switching 
between stable states of a periodically modulated oscilla- 
tor, we see that switching via tunneling can be observed 
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only where the relaxation rate is smaller than the tun- 
neling rate, that is the prefactor in the rate of quan- 
tum activation is very small. Such experiment requires 
preparing the system in one of the wells. Tunneling be- 
tween period two states of a strongly modulated strongly 
nonlinear system has been demonstrated in atomic optics 
[IE, m, For a weakly nonlinear oscillator, tunneling 
splitting between the lowest states of the Hamiltonian g 
for ^ ~ Q was found in Ref. [i^ . 

The effective Hamiltonian g{P, Q) is quartic in the 
momentum P. Therefore the semiclassical momentum 
P((5, g) as given by the equation Q) = g has 4 rather 
than 2 branches. This leads to new features of tunnel- 
ing compared to the standard picture for one-dimensional 
systems with a time- independent Hamiltonian quadratic 
in P. 

For concrcteness, we will consider tunneling from the 
left well, which is located at Q = Qio = —(1 -I- fJ,y^'^, P = 
0. The rate of tunneling from the bottom of the well of 
g{P, Q) in Fig. 2] is determined by the WKB wave func- 
tion with g = (7„iin = — (1 -I- At)^/4- This wave function is 
particularly simple for /i < 0. In the region Qiq < Q and 
not too close to Q/o it has the form 



iVz) =C(ap.g)^i/2gxp[z5o(Q)/A], 
rQ 

5o(g)= / P-{Q')dQ', 
Jqw 



where 



1/2- 



(41) 



(42) 



For /i < the wave function \tpi) monotonically decays 
with increasing Q. The exponent for interwell tunneling 
is S'tun/A, with 5'tun = Im So{-Qio) (we use that -Qio 
is the position of the right well), 



St, 



A more interesting situation arises in the case ^ > 



1/2 



< 



0. Here, inside the classically forbidden region — /i 
Q < y}/"^ decay of the wave function is accompanied by 
oscillations. We leave the analysis of this behavior for 
a separate paper. Here we only note that the tunneling 
exponent is still given by Eq. 

If the level splitting due to tunneling is small compared 
to their broadening due to relaxation, the tunneling prob- 
ability is quadratic in the tunneling amplitude. 



T^tun oc exp(-2S'tun/A). 



(44) 



The action 2S'tun as a function of /i is plotted in Fig. [TOl 
It is seen from this figure that the tunneling exponent 
exceeds the quantum activation exponent Ra/^ for all 
values of the control parameter ^. This indicates that, 
as mentioned above, it is exponentially more probable to 
switch between the classically stable states of the oscil- 
lator via activation than via tunneling from the ground 
intrawell state, for not too small relaxation rate. 




FIG. 10: (Color online) The scaled exponent 2Stun in the 
tunneling probability as a function of the parameter fi — 
2u}pSuj/F (solid line). Also shown for comparison is the quan- 
tum activation energy Ra for fi = and n ^ 0. 



In the same limit where the relaxation rate exceeds 
the tunneling rate, we can consider the effect of tunnel- 
ing from excited intrawell states of the Hamiltonian g. 
The analysis is similar to that for systems in thermal 
equilibrium [H, [3i 13 • Over the relaxation time there is 
formed a quasiequilibrium distribution over the states in- 
side the initially occupied well of g{P, Q). As before, we 
assume that relaxational broadening of the levels gn is 
small compared to the level spacing. The probability of 
tunneling from a given state n is determined by its occu- 
pation pn- The overall switching probability is given by 
a sum of tunneling probabilities from individual intrawell 
states 



■2S„/A 



Pn 



(45) 



where C„ is a prefactor that smoothly depends on ?i, and 



Sn = Im j P{Q, gn)dQ 



is the imaginary part of the action of a classical particle 
with Hamiltonian g{P, Q) and energy g„, which moves in 
complex time from the turning point P = in one well 
of g{P, Q) to the turning point in the other well. It fol- 
lows from the analysis in Appendix A and Eq. (|A3[) , that 

the duration of such motion is {t^^ -\- t^^)/2. Therefore 

c'5„/9g„ = -ImT^'V2. 

Taking into account that pn = exp(— i?„/A), we see 
that the derivative of the overall exponent in Eq. ([45]) 

over gn is Im —R! {gn). It was shown in Secs. lIVI aiidlVl 
that this derivative is always positive, cf. Fig. [51 There- 
fore the exponent monotonically increases (decreases in 
absolute value) with increasing g. This shows that, with 



13 



overwhelming probability, switching occurs via overbar- 
rier transitions, i.e., the switching mechanism is quantum 
activation. We emphasize that this result is independent 
of the bath temperature. 



IX. CONCLUSIONS 

In this paper we studied switching between the states 
of period two vibrations of a parametrically modulated 
nonlinear oscillator and the distribution over period-two 
quasienergy levels. We considered a semiclassical case 
where the wells of the scaled oscillator Hamiltonian in 
the rotating frame g{P,Q), Fig. [H contain many levels. 
The distance between the levels is small compared to 
hiop but much larger than the tunneling splitting. We 
assumed that the oscillator is underdamped, so that the 
interlevel distance exceeds their width. At the same time, 
of primary interest was the case where this width exceeds 
the tunneling splitting. 

We considered relaxation due to coupling to a thermal 
bath, which is linear in the oscillator coordinate, as well 
as dephasing from random noise that modulates the os- 
cillator frequency. The problem of the distribution over 
intrawell states was reduced to a balance equation. The 
coefficients in this equation were obtained explicitly in 
the WKB approximation, using the analytical properties 
of the solution of the classical equations of motion. The 
balance equation was then solved in the eikonal approx- 
imation. The eikonal solution was confirmed by a full 
numerical solution of the balance equation, which did 
not use the WKB approximation for transition matrix 
elements. 

We found that the distribution over period-two 
quasienergy has a form of the Boltzmann distribu- 
tion with effective temperature that depends on the 
quasienergy. This temperature remains nonzero even 
where the temperature of the thermal bath T ^ 0. It 
is determined by diffusion over quasienergy, which ac- 
companies relaxation and has quantum origin: it is due 
to the Floquet wave functions being combinations of the 
Fock wave functions of the oscillator. 

Unexpectedly, we found that the quasienergy distribu- 
tion for T = 0, where the system has detailed balance, 
is fragile. It differs significantly from the solution for 
T — * 0, for a large number of intrawell states. The T ~ 
solution is also destroyed by even small dephasing. 

The probability of switching between period two states 
Wsw is determined by the occupation of the states near 
the barrier top of the effective Hamiltonian g{P, Q). We 
calculated the effective quantum activation energy Ra 
which gives W^sw « exp(— _R^/A). Both Ra/^ and the 
exponent of the tunneling probability are proportional 
to the reciprocal scaled Planck constant A^^. However, 
for all parameter values and all bath temperatures, Ra/^ 
is smaller than the tunneling exponent. Therefore in the 
case where intrawell relaxation is faster than interwell 
tunneling, switching occurs via quantum activation. 



In the limit where fluctuations of the oscillator are 
classical, kT ^ hojp, we have Ra oc (kT)^^, and Wsw 
is described by the standard activation law. Down to 
small Planck numbers n > 0.1 the quantum activation 
energy Ra is reasonably well described by the classi- 
cal expression even for small kT /hwp provided kT is re- 
placed by ftwi?(2n4-l)/4, with n being the Planck number 
of the oscillator. The inapplicability of this description 
for small n indicates, however, that classical and quan- 
tum fluctuations do not simply add up. The replacement 
kT huipi'^n + l)/4 becomes exact for all fi close to 
the bifurcational value of the control parameter ^ = 
where the period-two states first emerge. In this range 
Ra scales with the distance to the bifurcation point as 
i?A oc - ^IbY- 

The results on switching rate are accessible to direct 
experimental studies in currently studied nano- and mi- 
crosystems, in particular in systems based on Joseplison 
junctions, including those used for highly sensitive quan- 
tum measurements. 

We gratefully acknowledge a discussion with M. De- 
voret. This work was partly supported by the NSF 
through grant No. ITR-0085922. 



APPENDIX A: CLASSICAL MOTION OF THE 
PARAMETRIC OSCILLATOR 

The Hamiltonian g{P^ Q) ^ is quartic in P and Q. 
This makes it possible to solve classical equations of mo- 
tion Trajectories with given g lie on the cross-section 
Q)= g oi the surface g(P, Q) in Fig. H We will be 
interested only in the intrawell trajectories, in which case 
g < 0. The trajectories in different wells are inversely 
symmetrical and can be obtained by the transformation 
Q ^ —Q, P ^ —P. Their time dependence can be ex- 
pressed in terms of the Jacobi elliptic functions [s^ . For 
trajectories in the right well in Fig. 21 where Q > 0, we 
have 



Q{r) 



P{r) 



23/2 |g| 1/2 jj^^/ 

+ K_cnr' 
Igj-'-Z'^snr' 
+ K-Cnr' 



(Al) 



Here 



K± = 1 + M±2|5 



1/2 



1/2 



23/2|.g|i/V. (A2) 



The parameter of the elliptic functions mj = mj{g) is 
given by Eq. (jl2p . For /i < the function mj{g) mono- 
tonically decreases with increasing g from uij = for 
9 = .9min = — + 1)^/4 to mj — > — oo for (7 — > 0. For 
/i > the function rnj{g) becomes nonmonotonic. It 
first increases from mj = with increasing g, but than 
decreases, goes trough mj = for 5 = —(1 — ^)^/4, and 
goes to —00 for g — > 0. The Jacobi functions sn(r'|mj), 
cn(r'|mj) and dn(r'|mj) for mj < are equal to 
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(1 — mj)~-^/^sd (f'|mj), cd(f'|mj), and nd (fl ffij ) with 
rhj = — mj/(l — mj) and f' = (1 — tuj^^'^t' [SSj . 

The double periodicity of the functions Q(t), P(r) dis- 
cussed in Sec. IIIBl is a consequence of the double period- 
icity of elliptic functions. The expressions for the periods 
T^^'^\ Eq. dm), follow from Eq. (|Ai|). The trajectories 
in the left well of g{P, Q) in Fig. 21 where Q < 0, can be 
written in the form 



r(2) 



(A3) 



Pi{t) = p{t 



-(1) 



r(2) 



/2 



This expression shows how to make a transition from 
one well to another by moving in complex time, which 
simplifies the analysis of oscillator tunneling. 

We note that the function cut' in the expressions 

for Q,P dXH) has periods t^^\ (t^^^ +r^^^) /2 [H as 

a function of r. It's period parallelogram is shown in 
Fig. [5l In this parallelogram cnr' takes on any value 
twice. Therefore both Q{t) and P{t) have two poles lo- 
cated at T = r*,r*». The values of r*,r,t, are given by 
the equation 



(2) 



- - ( ^(1) 
~ 2 VP 



(A4) 



p 



The positions of the poles are shown in Fig.[5]for mj > 
and mj < 0, respectively, with 0, — 27rr»/rp^'', = 
27rr**/Tp . For concreteness we choose 



Imr* < Imr^^^ /4. 



(A5) 



Using the relations between the Jacobi elliptic functions 
[38| we find that, near the pole at r = t*. 



(A6) 



P-^Q(x -(t-t,)-\ 

whereas P — iQ is not singular at t = t»». Eq. (jA6p was 
used to obtain the explicit form of the matrix element of 
the operator a = {2X)-^^'^{P - iQ) in Sec. II C. 

APPENDIX B: THE CLASSICAL LIMIT 



ber n ^ 1. The explicit form of the coefficients in 
Eq. (P7)) for R' follows from the general expression 
for the transition rates Wnn' , 



^■mWn+„in= ^ m|a_™(5„)p, (Bl) 

m 7n— — oo 

oo 

Wn+,nn = (2n + 1) X! '7l^|a-m(gn)P- 
m m— — oo 

The semiclassical matrix elements a_m are given by Eq. 
psp . They are Fourier components of the function a(<^; g) 
on the classical trajectory with given g = Qn- 



E 



Using 



the 

Jm(<p2-'t>l) 



completeness 



condition 



m— — oo 



2tt6{4>2 — 4>i) we can rewrite 



oo ^ „27r 

^ m\a-yn{9)? ^ ^ j d(l)a{(j)]g)d^a*{(l)\g), 

m=-oo "^O 
oo 

5] m2|a_™(g)p = (2^)-i 

(B2) 



m— — OO 



X / d(f>d^a{(t>;g)d^a*{(l);g), 
Jo 



with a(0;g) = {2Xr^/%P{T; g) - *Q(r;g)) and = 
Lu(g)T. The integrals over can be written as contour 
integrals over dP, dQ along the trajectories with given 
g. The contour integrals can be further simplified using 
the Stokes theorem. This gives 



(9)1' 



m— — OO 



2A7r 



(B3) 



m— — OO 



2A7r 



Ni9), 



In this Appendix we calculate the effective inverse tem- 
perature R'{g) in the limit of large oscillator Planck num- 



where the functions M{g) and N{g) are given by Eq. 

mi- 
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